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INVESTIGATION OF SONIC BOOM GENERATED BY THIN, 
NONLIFTING, RECTANGULAR WINGS 
Sanford S. Davis 
Ames Research Center 

SUMMARY 


A new theory is described for predicting sonic boom pressure signatures 
emitted by nonlifting rectangular wings. Comparisons are made with previous 
(Whitham) theory and with experimentally determined near field signatures. 
Although both theories agree with experiment for low aspect ratio wings , 
Whitham* s theory seriously overpredicts signature length for high aspect ratio 
wings. No experiments were conducted in the far field, but the two theories 
predict nearly identical results in this region. 


INTRODUCTION 


The standard methods used to predict the sonic boom are based on a theory 
developed by Whitham (refs. 1 and 2). Although the theory was originally 
devised for slender, axisymmetric configurations, Whitham extended his theory 
to include arbitrary configurations by using the methods of geometrical 
acoustics (ref. 3). He later realized that the procedures developed for 
slender, symmetrical bodies of revolution may be used intact for a nonsymmet- 
rical configuration if the so-called equivalent body of revolution is intro- 
duced. But if the equivalent body of revolution is not slender, Whitham' s 
theory as given in reference 2 is not strictly applicable. The extent to 
which the pressure signatures calculated for a nons lender equivalent body 
deviates from the experimentally determined signatures will be explored in 
this report. 

In order to predict sonic boom pressure signatures more accurately in the 
nearfields and midfields for nonslender configurations, a new theory has been 
developed. Although restricted to nonlifting, rectangular wings, this theory 
predicts substantial differences in pressure signatures when compared with the 
equivalent body theory. The purpose of this report is to evaluate signatures 
predicted by both theories and to compare them with experimentally determined 
signatures. 

The rectangular wing represents a very interesting model configuration 
because of the diversity of flows surrounding it. The wing tip Mach cones 
divide the flow into three distinct regions: (1) the two-dimensional region, 

(2) the tip- cone region, and (3) the interaction region. In the two- 
dimensional region, the flow is independent of the spanwise coordinate. The 
flow field is exactly the same as that generated by a two-dimensional wing 



with the same airfoil section as the rectangular wing. The tip-cone region 
contains a three-dimensional flow, a small portion of which has conical symme- 
try. In the interaction region the flow is fully three-dimensional. There- 
fore, a theory that predicts the shock wave generated by a rectangular wing 
includes elements of the theories of conical flows, two- and three-dimensional 
flows, and the transition flows that separate these regions. It will be shown 
that the new theory, called the uniform theory, more correctly describes the 
flow in each of these regions. 

The flow geometry and the uniform theory are reviewed briefly in the 
following sections. The theory is developed for a parabolic arc airfoil sec- 
tion, but the uniform theory can be applied to any wing with a biconvex 
airfoil section. 


LINEAR FLOW FIELD GEOMETRY 


Consider the wing to be fixed in an (x,y,z) Cartesian coordinate system. 
The mean plane of the wing occupies part of the x-y plane in the region 
x > 0, y > 0. A supersonic stream flows in the direction of increasing x at 
Mach number M. Since the flow is symmetric with respect to the plane z = 0, 
only the region z > 0 will be considered. 


The thickness distribution of the wing is denoted by t(x) and is given 
by the following formula: 


t(x) = 0 

all y 

X 

< 0 

t(x) = 2m (x - x 2 /2d) 

0 < y < l \ 

0 

< x < 2d 

t(x) = 0 

y < 0 , y > l 1 



t(x) =0 

all y 

X 

> 2d 


( 1 ) 


In equation (1) £ represents the span of the wing, 2d is its chord, and 
m is the slope of the airfoil section at the leading edge. Figure 1 is a 
plan view of the wing. 


A fundamental property of supersonic flows is the separation of the flow 
field into disturbed and undisturbed regions. In linearized theory, the 
boundary between these two regions is a characteristic surface (i.e., the wave 
front). For the rectangular wing, the linearized wave front consists of the 
following surface in z > 0. 


1. The planar surface 


where 


x - gz = 0 0 < y < l 

3 2 _ _ 


2 


1 



2. A portion of the tip cone at (0, 0, 0) 


x - 3 2 i/y 2 + z 2 = 0 y < 0 

3. A portion of the tip cone at (0, Z, 0) 

x - $/(y - Z -) 2 + z 2 = 0 y > Z- 

If the tip cones, given by 2 and 3 above, are extended to y > 0, y < Z, 
three distinct disturbed regions are defined. These regions are denoted in 
figure 1 by the numerals I, II, and III. In region I the flow is two- 
dimensional, the same as for a wing of infinite span. In region II, inside 
one or the other of the tip cones, the flow is dependent on only one tip cone. 
Regions I and II alone represent the flow field for a semi-infinite wing. 
Region III is the interaction region. In this region the flow field is 
influenced by both tip cones. Region III represents the typical flow field 
pattern at large distances from a finite wing. 

The general appearance of these regions, at various distances from the 
wing, is also shown in figure 1. In each of the sketches labeled nearfield, 
midfield, and farfield, the wave front is shown as it would appear when inter- 
sected by a plane at a constant value of x. In the nearfield sketch, which 
represents a distance downstream of one or two spans, regions I, II, and III 
are about the same size. In the midfield, at a few spans downstream, 
regions I and II are reduced at the expense of region III. Many spans down- 
stream, the characteristic farfield pattern emerges. Region I becomes very 
small, region II becomes a crescent shaped area with a thickness the order of 
the span, while region III dominates the flow. 


DESCRIPTION OF THE UNIFORM THEORY 


The basis of the uniform theory is Whitham's hypothesis (as developed 
and applied in refs. 2 and 3), which states that linear theory gives a correct 
first approximation to variations of the physical quantities along the lin- 
earized (free-stream) characteristics, but the location of these linearized 
characteristics is in error. In order to apply Whitham's hypothesis to the 
rectangular wing, a better approximation to linearized theory must be obtained 
than the one afforded by the equivalent body method. Unfortunately, a better 
approximation cannot be expressed in the form of the Whitham F- function. Once 
this form of solution is abandoned, Whitham's simple and elegant results can- 
not be duplicated. But there is no reason why an improved form of the linear 
solution could not be used as the basis of a characteristic correction, and 
this in fact is the crux of the uniform theory. This theory will be 
described briefly in the following paragraphs. 

An expression for the perturbation velocity potential due to the presence 
of a nonlifting rectangular wing in an otherwise undisturbed stream of Mach 
number M is 


3 


m(l - x 1 /d)dx 1 dy 1 


( 2 ) 


4> 


V 

ir 



»/(x - xi) 2 - 3 2 (y - yi) 2 - g 2 z 2 


In this expression m(l - xi/d) represents the streamwise slope distribution 
of the parabolic airfoil, V is the free-stream velocity, and A is the area 
in the plane of the wing bounded by the Mach forecone from the field point 
(x,y,z), and the edges of the wing. 

In the region z > 0 , 0 < y < l a new set of coordinates based on the 
family of planar characteristics will be introduced. Let 

x = + Sri 

X = Pi 
z = rj 


If the location of the point Pi in figure 2 is expressed in terms of 
(£ x , ri, p x ), the potential can be written as 


$ = 


V 

ir 



m(l - x i /d)dx 1 dy 1 
______ 


(3) 


= (Ci - xi) 2 + 2gri(5i - xi) - S 2 (pi - yi) 2 


In the region near the wave front, £i is always much smaller than ri (see 
fig. 2).. Since only the region near the front is of interest, equation (3) 
may be simplified by neglecting the first term in the expression for Xi . If 
this quadratic term in (£i - Xi) is neglected, the resulting expression for $ 
is valid for small values of ?i/ r i* The simplified form of equation (3) is 


$ 


<f> = 


; /I 


m(l - xi/d)dx 1 dy 1 


A v4eriUi - xi) - e 2 (pi - yj 2 


(4) 


where A now represents the area in the plane of the wing bounded by 


2gri(C 1 - x x ) - g 2 (p 1 - y x ) 2 = 0 

and the wing edges. Note that in this approximation the intersection of the 
plane of the wing and the Mach forecone is represented by an approximating 
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parabola instead of a hyperbola. (The F- function approach, by neglecting the 
quadratic term B 2 (p 1 - yj) 2 , would represent this approximation by a straight 
line parallel to the leading edge.) Equation (4) represents the uniformly 
valid approximation to the linear perturbation potential in the region z > 0, 
0 < y < l. 

Another set of coordinates is chosen in the region z > 0, y < 0. (The 
region z > 0, y > l will not be examined since the plane y = 1/2 
represents a plane of symmetry.) 

x = 5 2 + 3r 2 
y = r 2 cos p 2 
z = r 2 sin p 2 

If the location of the point P 2 in figure 2 is expressed in terms of (£ 2 , 
r 2 , p 2 ) , the potential in this coordinate system is 


V |*r -d " x 1 /d)dx 1 dy 1 


X 2 = (S 2 - Xi ) 2 + 2Br 2 (£ 2 - x x + 0y x cos p 2 ) - B 2 y x 2 

Equation (5) can be approximated in the region where £ 2 / r 2 sma1 ^ if the 
term (£ 2 - x l) 2 i n ^2 i s neglected. This approximate potential is written 


$ 



m(l - xi/d)dxidyi 

,/ 23r 2 (5 2 - X! + Byj cos p 2 ) - B 2 y x 2 


where A represents the area in the plane of the wing bounded by 


( 6 ) 


2|3r 2 (? 2 - xi + Byi cos p 2 ) - 3 2 yj 2 = 0 


and the wing edges. As mentioned above, this approximation results in a 
parabolic forecone intersection. (It does not give as crude an approximation 
as the F-function in which X 2 is linearized by neglecting 3 2 yi . ) Equa- 
tion (6) represents the uniformly valid approximation to the potential in the 
region z > 0, y < 0. It can be shown that this form of <j> is a uniformly 
valid asymptotic representation to $ with respect to the parameter p 2 as 
? 2/ r 2 0 . 
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The choice of new coordinates in equations (4) and (6) is not arbitrary, 
but is closely related to the rays and wave fronts of geometrical acoustics. 
In this context, the x-axis plays the role of time and the y-z plane repre- 
sents the space; rj and r 2 represent distance along a "ray"; pi and p 2 are 
ray labelling parameters; and ^ and £ 2 are "phase variables." Flow past a 
rectangular wing is analogous in many respects to diffraction of a pulse by a 
plane wall, and "shadow regions" or "reflected regions" may be defined in the 
y-z plane (e.g. , region I is a reflected zone, region II is a shadow region 
and the shadow boundary is the ray pj = 0 or p 2 = tt/ 2). This analogy has 
been very useful in determining the expressions (4) and (6) . 

The approximations (4) and (6) may be integrated directly, resulting in 
the following expressions for the axial perturbation velocity: 

1. In the region z>0, 0 < y < l: 


V - it D4)<" - •' *D- l&4- S( Vdrf AZ ] | 


/ 2£ 2 r^ 

A ° " Re v & Pi 2 ’ l . 


1/2 


A i = Re 


r ^ -ii 

L B(Z - p x ) 2 J 


1/2 


(7) 


2. In the region z > 0, y < 0: 


u -m \L h er 2 cos2 -i„ * -1 n , 

V = ¥6 - T 2d ) (tan B o “ tan 

gr 2 cos 2 0B O 3(r, cos 0 + Z) 2 1 

H\ 


Bn - Rfi 


B i = Re 


2d 

25, 


2r 2 d 


, 1/2 


3r 2 cos 2 0/ 

(2£ 2 r 2 /B) ~ 2t 2 1 cos 0 - l 2 
(r 2 cos 0 + V) 2 


“I 1/2 


( 8 ) 


where 0 = ir - p 2 . 



Note that equations (7) include a region where both A 0 and A i vanish. 
The perturbation velocity in this region is the same that would exist for a 
two-dimensional flow in the x-z plane. As A 0 and A£ become real, the 
tip cones modify this two-dimensional region. At very large values of r^ , if 
?! 0, the solution reduces to the form that would be obtained by equivalent 

body theory. Likewise the form given by equations (8) reduces to the equiva- 
lent body solution as r 2 -*■ °°. (See the appendix for the equivalent body 
corresponding to this wing.) 

Returning to Whitham' s hypothesis, the linearized solution given by 
equations (7) and (8) must be corrected to account for cumulative nonlinear 
effects. The correction procedure is to redefine the lines = x - 8rj and 
S 2 = x ■ 3r 2 to be the next approximations to the exact characteristic lines. 
The slope of the next higher-order characteristic lines in the planes 
p^ = constant is 


dx _ „ . (y + 1)M 4 u (5i> r i*Pi) w9 ["„ u (5i>*i»Pj.) Ur i^i* r i*Pi) 

d?T ' 3 + 2 $ 2 v M [_ e v + v 


+ terms of order 





(9) 


where i = 1 for z>0, 0 < y < l and i = 2 for z > 0, y < 0. In equa- 
tion (9) Uri/V represents the velocity in the direction of increasing r^ 
and y = Cp/c v is the ratio of specific heats. Since only a first-order 
characteristic correction is required, only the linear terms in the exact 
expression for the characteristics have been included in equation (9). In 
addition, the term 8(u/V) + (ur-j/V) may be neglected because the curvature of 
the characteristic line is very small. (For planar families of characteristics 
8(u/V) + (u^/V) is identically zero.) 


Equation (9) may be integrated, and, in analogy with the procedure 
developed by Whitham, the corrected characteristics, may be expressed as 


x - 8ri - G(? i ,ri,p i ) + q 


GUj 


i» r i>P i ) = - 2 b 2 )M4 1 V (5i,t,pi)dt 


( 10 ) 


The function G(?^,rj_,p i ) is given explicitly as: 
1. In the region z>0, 0<y<Z 
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GCSx.rj.Pj) = 


(Y 


a* 5 [(■ - $ (-. - » 


8p 2 A q 

tan -1 A 0 + — *tan" 1 A^ 


2 . 


3CZ--P i) 2 A^\ Bpj 2 
+ 25i /' d 

In the region z > 0, y < 0 


B(Z - p,) 2 

j j— (Ao - tan" 1 A 0 ) 3 (A^ - tan" 1 A^j 




-ro tan 1 B 7 — ^ 

^ ^ cos 9 


-fj^2 tan J B 0 + r 2 B Q - 
Tr, cos 0 + Z 2Co 


2€, 

5 — tan 1 B n 

B cos 2 0 0 


B Z 


B cos 2 0 


+tan 


r 2 cos 0 + Z 
BZ ^ B 2 


[ B cos 0(r 2 cos 0 + Z) 

2? 2 - BZ cos 0 B zJj ~ d 

BZ 2 -1 H r 2 cos 6 + «»zl 6r ? 2 cos2 6 , 

+ — : — tan 1 ^ 1 - 


Br 2 cos 2 0 Br 2 cos 2 0 „ pi „ 

2 n 2 R 3 2 


4d 


Bq + 


4d 

Br, 2 cos 2 0 


tan" 1 B r 


i 


4d 


4d 


Bq 4 tan" 1 B -1 


Br 2 2 cos 2 0 B 2 Z 2 cos 0(r 2 cos 9 + Z) 

4d tan B Z- + 8d(5 2 - BZ cos 0) B Z 

£ 2 (2£ 2 - gZ cos 0) (r 2 cos 0 + Z) 

4d cos 0(£ 2 " BZ cos 0) B ^ 


B 2 cos 0 (r 2 cos 0 + Z) 
8d(£ - BZ cos 0) 


B z 3 


Br . 2 cos 2 0 
+ 4 ^ B Q tan 


1 f' 


B cos 0(r 2 cos 0 
2£ 2 - BZ cos 0 


^, 1 ) 


Equations (7) , with £4 given by equation (10) represent the uniformly valid 
first-order expression for the velocity field in the region z > 0 , 0 < y < Z. 
Equations ( 8 ), with £ 2 given by equation (10), represent the analogous 
result for z > 0 , y < 0 . 

8 



The solutions represented by equations (7) and (10, i=l) or equations (8) 
and (10, i=2) are multivalued, and thus do not represent a physically plausi- 
ble situation. In order to eliminate the multivalues, a shock wave must be 
introduced into the flow field. The shock wave separates the regions of dis- 
turbed and undisturbed flow, and is located by the rule that the direction of 
the shock wave bisects the angle between the characteristics which meet at a 
point. 

In the undisturbed flow, the slope of the characteristic is the 
free-stream wave slope 8. In the disturbed flow, the characteristic direc- 
tion to be used is associated with the family that carries local discontinu- 
ities in the perturbation quantities. In the case of the shock produced by a 
rectangular wing, more than one such family exists. For example, figure 3 
shows the distorted characteristics in the plane of symmetry of the wing. (In 
the plane of symmetry (pj = 1/2) both tip cones have the same trace.) In the 
region ahead of the trace of the distorted tip cone characteristic, the family 
of distorted planar characteristics is used to calculate the shock direction. 
But within the influence domain of the tip cone, the family associated with 
the tip cone trace is used to determine the shock-wave slope. 

Each family of characteristics may be expressed as 


= ? i - f CV 


where £ 0 identifies the family to be used in the bisection rule, and f(r^) 
is a function appropriate to each family. In figure 3, where i = 1, the two 
families used are (1) £ 0 = £., f(ri) = 0 and (2) = £, - 8(£/2) 2 /2r,, 

f( ri ) = B(^/2) 2 /2r 1 . 

In terms of the parameter £ 0 , the corrected characteristics can be 
expressed as 


x = $ri - G[£ 0 + f (r-jJ ,r^] + K 0 + f(r A ) 


( 11 ) 


The bisection rule, when applied to (11), becomes 


Slope of shock = J [(3) * (a - 


C1 - . . 0 1 3G 1 3f 

Slope of shock = g - - — + - — 


If the position of the shock wave is given by 

x sh = &*i " H ( r i) 


( 12 ) 


(13) 
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then the slope of the shock dx^/dr^ is 


dx sh 

dri 


= 8 - 


dH(r.) 

dri 


or 


dH _ _1 3G_ _ 1 3f 
dri ~ 2 3r^ " 2 3r^ 


(14) 


(15) 


Note that all partial derivatives in equation (15) are to be taken with £ 0 
held constant. 

From equations (11) and (13), H may be expressed as 


H = G - S o - f (16) 

Equations (15) and (16) represent a first-order, nonlinear, ordinary differ- 
ential equation for f; as a function of r-^. The functional form of G is 
so complicated that a numerical approach must be used to solve this equation. 

In order to solve equations (15) and (16) a modification of the familiar 
Euler-Cauchy forward integration procedure was used. In terms of small, 
finite increments the differential system can be written as: 


AH = _1 3G_ _ 1 3f_ 
Ari 2 3r^ ~ 2 3r^ 


(17a) 


AH = AG - A£ 0 - Af 


(17b) 


If starting values £ 00 and r^ Q are chosen, and an increment Ar^ is given, 
equation (17a) shows that H changes by an amount 


AH 


(l 3G_ 
\2 3ri 


1 M_ 

2 3ri 


)l 


Ar-; 


^oo ,r io 


(18a) 


When r^ = r^ 0 + Ari,£ Q must change to £ 00 + A£ q . Equation (17b) reflects 
these increments as a change in AH by an amount 


AH = G(£ 00 + A? 0 ,r i0 + Ar ± ) - GU oQ ,r. o ) - A? q - f(r io + Ar f ) + f(r io ) 


(18b) 
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Since H is the same in both equations (18a) and (18b) , 


1 3G_ 

2 3ri 


1 3f_ 

2 3ri 


^oo ,r io 


Ar i = G (5 00 + A So> r io + Ar i) " G (5 00 ,r i0 ) 

- f ( r io + A^) + f ( r i 0 5 


A ?o 


( 19 ) 


Equation (19) is a nonlinear algebraic equation for determining the unknown 
increment A£ 0 . Once A5 0 is found, a new starting point 


5o = 5 00 + A 5 0 ’ r i = r iQ + Ar ± 


is found. The procedure is continued in this manner until the shock wave 
intersects a new family of characteristics. At this point £ 0 is redefined, 
the starting value of r 0 is obtained from the preceding shock wave calcula- 
tion, and the numerical integration can be restarted. 

A slight difficulty in the procedure appears in the region z > 0, y < 0. 
For the starting values, 5 0 = 0, = 0 when substituted into the governing 

formulas, gives a result which is the indeterminate from 0/0. This problem 
is easily surmounted by matching to a conical flow solution near T 2 = 0. 

Once the uniformly valid velocity field (with shock waves) has been 
obtained, the pressure rise at any point in the flow is calculated from the 
expression 


= - v M 2 " 

p Y V 

r CO 

TEST PROCEDURES AND RESULTS 


The general arrangement of the wing model is shown in figure 4. The 
symmetric airfoil section is a portion of a parabolic arc with a thickness 
distribution given by t = 0.250 (x - 0.197 x 2 ) cm where t is the total 
thickness of the wing and x is a streamwise variable. Wings were fabricated 
with spans of 2.54 cm (Z/ d=l), 10.16 cm (Z/d = 4) and 15.24 cm (Z/d = 6). 

The, experiments were conducted in the Ames 2- by 2-Foot Transonic Wind 
Tunnel at a Mach number of 1.40, a total pressure of 67.5 kN/m 2 and a total 
temperature 296° K. The arrangement of the model in the wind tunnel is shown 
in figure 5. The wing is supported in the test section by a tapered, conical 
sting faired into the wing’s axis of symmetry. In order to minimize its 
effect on the front shock wave, the sting was fastened to the aft portion of 
the wing only. The 58.5 cm long sting assembly was secured to a linear 
actuator. The actuator had a continuously variable 15.24 cm streamwise motion. 
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and could be preset to any streamwise station prior to a test run. During a 
run, the x (streamwise) position of the wing was monitored by a calibrated 
potentiometer that was coupled to the 15.24 cm motion. 

The shock-wave signature was measured with a slender 1° semiangle probe. 
The probe was offset 10.2 cm from the test section wall in order to exclude 
the effect of the wall's boundary layer from the measurements. On a plane 
inclined at 45° to the flow direction, four holes of 0.1024 cm diameter were 
drilled 7.62 cm from the probe tip. The pressures from this probe were 
measured relative to the free-stream static pressure by means of a calibrated 
pressure transducer. 

Before each run, the wing was set in the test section so that the shock 
wave would be slightly behind the probe orifice. During the run the model was 
moved upstream in discrete steps by using the 15.24 cm motion of the actuator. 
At each step, the overpressure and distance moved was recorded. In this 
manner the entire shock-wave signature was obtained. 

The results of the test program are summarized in figures 6 through 8. 

In each case the signatures are shown as normalized overpressures as a func- 
tion of a horizontal scale measured from the position of the linearized lead- 
ing Mach line. In all the figures, linear dimensions have been scaled with 
respect to the distance from the leading edge to the point of maximum 
thickness (i.e., with respect to 2.54 cm). 

Figure 6 shows the shock-wave signatures for each of the three wings 
tested at varying distances below the wing (z direction) . The wings having 
values of 7/d of 4 and 6 show strong two-dimensional shock waves in the near 
field, while the 7/d = 1 wing shows a relatively weaker three-dimensional 
type of shock wave signature at all distances. The relative change in peak 
shock overpressure with distance is greatest for the wings with 7/d of 4 and 
6. This effect is attributed to the sharp pressure relief associated with the 
wing tip cones. This behavior is indicated in figure 7(a) which shows the 
shock signature for the wing 7/d = 4 and z/d = 2.2, but at different span- 
wise stations. The shock strength below the wing tip (span station y/d = 0.0) 
is about half the value at the centerline (span station y/d = 2.0) because 
the centerline station is in a two-dimensional flow field (as far as the shock 
strength is concerned), while the tip station is located within the dependence 
domain of the tip cone at y/d = 0.0. This ratio of 2:1 in peak overpressure 
is also predicted by the linearized theory, but (of course) the magnitudes 
obtained from linear theory are incorrect. At larger distances, the effect of 
spanwise variation is reduced. This lateral smoothing of the pressure field 
is due to the net cancelling effect of both wing tip cones on the two- 
dimensional field. This effect is shown in figure 7(b). 

The influence of roll angle on the measured pressure signatures is shown 
in figure 8 for each of the three wings. In each case, the roll angle is 
defined by a rotation of the plane of the wing relative to the original 
horizontal plane about the axis of symmetry of the wing. (This axis is also 
the axis of the sting.) Distance from the axis of symmetry is denoted by R/d, 
where 
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f7(i-ji) 2 '(t) 2 

For the 1 / d = 1 wing, the signature is almost axisymmetric at R/d = 16.0 
because at this distance the dominant effect on the sonic boom for this low 
aspect ratio wing is an axial source distribution. For the other wings large 
differences in meridional shock strengths and locations are apparent. 


COMPARISON OF MEASUREMENTS AND THEORY 


Predicted signatures of the uniform theory and equivalent body of 
revolution theory are compared with experimental signatures in figures 9 to 11. 
Details of the equivalent body theory for this configuration are described in 
the appendix. The signatures in the plane of symmetry below the wing for 
various length ratios, 1/ d, are compared in figure 9. These signatures indi- 
cate that for short wings (i.e., l/d = 1, fig. 9(a)) both the uniform theory 
and the equivalent body theory predict the experimental results quite well at 
all the distance ratios shown. However, the signatures corresponding to 
longer wings (i.e., figs. 19(b) and 9(c)) indicate that while the equivalent 
body theory gives reasonable values for the peak overpressure, it substan- 
tially overestimates the signature length. On the other hand, with one excep- 
tion, the uniform theory gives a better prediction for both the peak pressure 
and signature length for all values of l/d. The case z/d = 4 in fig^- 
ure 9(b) is an exception because the tip-cone interaction has not been pre- 
dicted exactly. This local defect has little effect on the overpressure at 
greater distances. The experimental signatures for z/d =8.0 and 16.0 
(fig. 9(c)) are probably less accurate than the others because the larger wing 
vibrated quite visibly during the test. These vibrations would tend to dis- 
tort the signatures at the larger distance ratios. An indication of this 
distortion is the rapid expansion behind the shock wave for the case 
z/d = 16.0 as compared to that for z/d = 8.0. All previous experimental and 
theoretical results indicate that the slope of the curve in the expansion 
region decreases, rather than increases, in a regularly decaying shock wave 
signature. 

In the midfield and farfield regions no experimental results are 
available, but a comparison of the two theories in their prediction of the 
peak overpressures is shown in figure 10. For l/d = 1, the agreement is 
good at all distances. For l/d = 4 and l/d = 6, the two theories lose mutual 
agreement at short distances. The nearfield experimental peak overpressures 
are also shown in figure 10, and these show that except for the one point 
corresponding to the l/d = 6 wing, good agreement between the uniform theory 
and experiment is achieved. 

An important property of the decay curves shown in figure 10 is the 
behavior of the uniform theory curve for the higher span wings . These curves 
approach the equivalent body theory from below and show that some benefit may 
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be obtained by using this "non F- function” behavior in the midfield region. 
Further work on this point is necessary before any conclusions concerning the 
benefits of an unswept wing may be reached. 

So far only the sonic boom overpressures under the wing have been 
compared. Figure 11 shows how the overpressures compare in the region to the 
side of the wing at z/d = 16.0. At this angle and distance, the uniform 
theory and the equivalent body theory predict similar signature shapes with 
reasonable agreement between theory and experiment. 


CONCLUSIONS 


A new theory for predicting sonic boom pressure signatures produced by 
nonlifting rectangular wings has been described and compared with the equiva- 
lent body of revolution theory and with a series of experimentally determined 
nearfield signatures. For low aspect ratio wings (i.e., less than one-half), 
both theories are in excellent agreement with experiment. For higher aspect 
ratio wings, with signatures taken in the vertical plane of symmetry, the 
equivalent body theory seriously overpredicts signature lengths. However, it 
does predict peak overpressure with accuracy comparable to the new theory. 

With signatures taken in the horizontal plane of the wing, the two theories 
are in excellent mutual agreement, and collectively are in reasonable agree- 
ment with experiment. Although no experiments were conducted in the farfield, 
the two theories give essentially identical predictions in this region. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Sept. 2, 1971 
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APPENDIX 


EQUIVALENT BODY METHOD 


The equivalent body for a rectangular wing is different for each 
meridional angle 9. The proper form of the equivalent body can be obtained 
from equation (6) for the perturbation velocity potential. If the radical in 
equation (6) is linearized, the potential may be expressed as 


♦ = - 



m(l - x 1 /d)dx 1 dyj 
/ 2 gr 2 /?2 ~ X 1 ~ ftyi cos 6 


(Al) 


where A represents the area in the plane of the wing bounded by 

- xi - Byj cos 9=0 and the wing boundary. An oblique coordinate system 


s l = y l 

ti = xi + By 1 cos 0 


CA2) 


is introduced into equation (Al) . In terms of (si,ti) as variables of 
integration, the potential becomes 



(A3) 


where the integral in the square bracket is taken over all values of s^ on 
the wing. 


From reference 2, the approximate perturbation velocity potential for a 
slender body of revolution may be expressed as 


<P = ~ 




(A4) 


where S(tj) is the cross-sectional area of the body at station tj, and 
primes indicate differentiation with respect to tj. A comparison of equa- 
tions (A3) and (A4) shows that the disturbance field produced by the rectan- 
gular wing is the same as would exist for a slender body of revolution with a 
cross-sectional area distribution of 
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\ S'(ti) = I d Sl m[l - ^ + cos e] CA5) 

•'wing 

Equation (A5) is valid for all values of ti less than 2d (i.e., excluding 
the influence domain of the trailing edge) and S' (tx) has the following 
representations in 0 < tx < 2d. 

fti/e cos 0 r t /gs x 

2 J Q " t + \~r) cos e J 

2 f o ds i m f 1 - -r + (^r) cos e ] 

Equation (A6) can be integrated to give an explicit formula for the 
cross-sectional area distribution of the equivalent slender body of revolution. 





The radius distribution is given by R(tx) = (l/ir)/S(tx) . 

Figure 12 shows a sampling of typical equivalent bodies of revolution 
for the wing l/d =4. At 0=0° and 0 = 45° the equivalent bodies are 
slender, but at 0 = 90° the equivalent body degenerates to a blunt forebody. 
Once these equivalent bodies have been found, the sonic boom corresponding to 
each of them may be found in a straightforward manner by following the 
procedures given in reference 2. 
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Figure 3.— Shock wave and distorted characteristics in the plane P! = (1/2) X . 










(a) l/d = 1 


Figure 6.— Effect of distance on the shock-wave signature measured in the plane of symmetry, M = 1.4. 
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Figure 9.— Comparison of uniform theory, equivalent body theory, and experiment for the overpressure signature in the plane of 

symmetry, M = 1.4. 
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(a) Z/d = 1 

Figure 1 0.— Decay of peak pressure rise with distance below wing. 
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Figure 10.- Continued. 
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(c) 1 / d - 6 


Figure 10.— Concluded. 
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Figure 11.— Overpressure signatures measured at R/d = 16.0 for a roll angle about the axis of symmetry of 90°, M = 1.4 
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